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ABSTRACT 

Motivated by the current constraints on the epoch of reionisation from recent cosmic microwave background observations, ionising 
background measurements of star-forming galaxies, and low redshifts line-of-sight probes, we propose a new data-motivated parame- 
terisation of the history of the average ionisation fraction. This parameterisation describes a flexible redshift-asymmetric reionisation 
process in two regimes that is capable of fitting all the current constraints. 
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1. Introduction 

After recombination at z 1090, the universe was essen¬ 
tially neutral. Observations of the Gunn-Peterson (GP) effect 
(Gunn & Peterson 1965) in quasar spectra (e.g. Fan et al. 2006) 
indicate that intergalactic gas has become almost fully reionised 
by redshift z ~ 6. This transition, called cosmic reionisation, is 
one of the most important events in cosmic structure formation 
and in cosmology. It is also of particular importance since it may 
be related to many fundamental questions, such as the presence 
of annihilating or decaying dark matter particles, properties of 
the first galaxies, physics of (mini-)quasars, and the formation 
of very metal-poor stars. The reionisation process starts with 
the formation of the first generation of early star-forming galax¬ 
ies and quasars that emit ultraviolet radiation that reionises the 
neutral regions around them. This is the so-called patchy reion¬ 
isation stage (Aghanim et al. 1996). After a sufficient number 
of ionising sources have formed and ionised regions have over¬ 
lapped, the ionised portion of the gas in the Universe rapidly in¬ 
creases until hydrogen becomes fully ionised. This period, dur¬ 
ing which the cosmic gas went from neutral to ionised, is known 
as the Epoch of Reionisation (EoR). The transition from neutral 
to ionised diffuse gas is still not observed directly but is con¬ 
strained by the absorption spectra of very distant quasars and 
gamma ray bursts, which reveal neutral hydrogen in intergalac¬ 
tic clouds. These GP effect observations show that the diffuse 
gas in the Universe is mostly ionised up to a redshift of order 6 
(Fan et al. 2006). 

The reionisation process describes the balance between the 
recombination of free electrons with protons to form neutral hy¬ 
drogen and the ionisation of hydrogen atoms by photons with 
energies hv > 13.6 eV. The first empirical, analytic, and nu¬ 
merical models of the reionisation process (e.g. Aghanim et al. 
1996; Gruzinov&Hu 1998; Madau et al. 1999; Gnedin 2000; 
Ciardi et al. 2003) highlighted the basic physics that give rise 
to the ionised intergalactic medium (IGM) at late times and 
provided predictions of the effects on the cosmic microwave 
background (CMB). The reionisation history is conveniently ex¬ 
pressed in terms of the filling factor of ionised hydrogen Qhii- 
One of the relevant physical and the most commonly used quan¬ 
tities for characterising reionisation is the Thomson scattering 


optical depth to the CMB r = ariecrj dn, where rig is the 
number density of free electrons at a conformal time rj, cr^ is the 
Thomson scattering cross-section, a the scale factor, and t/q the 
conformal time today. 

Reionisation leaves imprints on the CMB power spectra, 
both in polarisation and in intensity, and through the kinetic 
Sunyaev-Zeldovich (kSZ) effect, due to the re-scattering of pho¬ 
tons off free electrons (see Aghanim et al. 2008, and references 
therein). The temperature power spectrum is damped at scales 
larger than the horizon. This effect is degenerate with the am¬ 
plitude As and slope Us of the initial power spectrum. The 
CMB temperature thus constrains the amplitude The 

large-scale polarisation spectra present a reionisation “bump” 
that WMAP first measured in TE in 2003 (Spergel et al. 2003). 
In fact, CMB observations allow us to shed light only on the 
Thomson optical depth r. In addition to the CMB, other astro- 
physical probes constrain reionisation. 

Direct measurements of the ionisation state of the Universe 
through Lyof emission or quasars, between redshift 5 and 
redshift 10, have improved in quality and statistics (e.g. 
Fan et al. 2006; Faisst et al. 2014). Ionising background esti¬ 
mates from star-forming galaxies improved as well in par¬ 
allel to the depth of observations (e.g. Ishigaki et al. 2015; 
Robertson et al. 2015; Bouwens et al. 2015, hereafter ISHI15, 
ROB 15, BOU15, respectively). New probes are becoming avail¬ 
able and increasingly useful, such as the 21cm brightness fluc¬ 
tuations (e.g. Bowman & Rogers 2010) or gamma-ray bursts 
(e.g. Chornock et al. 2014). While it was hard to reconcile the 
ionising power of galaxies with the high value of r derived 
from WMAP (Spergel et al. 2003), the latest results from Planck 
(Planck Collaboration XIII 2015, hereafter PCP15) give an al¬ 
most coherent picture of the EoR where star-forming galaxies 
play a dominant role in ionising the Universe. 

Large numbers of detailed models of the reionisation his¬ 
tory have been constructed to describe the transition from a neu¬ 
tral to an ionised state and to reproduce the EoR observables 
(Pritchard et al. 2010; Pandolfi et al. 2011; Mitraetal. 2011, 
2012, 2015). These physically motivated models need to make 
assumptions on various complex processes (e.g. radiative trans¬ 
fer, star formation and evolution, chemical enrichement) cap- 
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tured in effective parameters, such as escape fraction, clumping 
factor, and star formation efficiency. As succesful as they are, 
those models cannot be used easily to explore the global history 
of the reionisation transition rapidly and without any prior on the 
detailed physics. 

In this letter, therefore, we provide a new parameterisation 
of the mean ionised fraction of the universe as a function of 
redshift. We exploit in particular the recent results based on the 
measurements of UV and IR flux of early galaxies in order to 
derive a simple data-motivated, fast, and economic expression 
of the ionised fraction. It agrees with the current observational 
constraints on EoR better than the standard parametisation used 
in CMB analyses. 

2. Observational status and current constraints 
on reionisation 

The reionisation process is intimately linked to the formation 
and evolution of the first emitting sources. At redshifts as high 
as 15 to 30, low-mass halos of matter detach from the general ex¬ 
pansion, collapse, and virialise. If they reach high enough tem¬ 
peratures to excite the lowest electronic levels of hydrogen, they 
cool and eventually form Population III stars and dwarf atomic- 
cooling galaxies. These emitting sources, thought to have ini¬ 
tiated the reionisation process, have different various properties 
(masses, soft/hard spectra, UV escape fractions, etc.). The reion¬ 
isation process is thus quite complex since it implies multiple 
sources of UV photons that combine their relative efficiencies at 
different epochs and scales. 

The abundance of quasars declines rapidly beyond red- 
shift z ~ 6. Quasars therefore cannot be significant contrib¬ 
utors to reionisation at high redshifts (e.g. Willott et al. 2010; 
Fontanot et al. 2012), but they are key players at lower red¬ 
shifts since they ensure the second reionisation of Helium 
(McQuinn et al. 2009). Faint star-forming galaxies at redshifts 
z > 6 have therefore been postulated as the most likely 
sources of reionisation, and their time-dependent abundance 
and spectral properties are crucial ingredients for understanding 
how intergalactic hydrogen became reionised (for reviews, see 
Fan et al. 2006; Robertson et al. 2010; Foeb & Furlanetto 2013). 
The high-redshift galaxy UV luminosity function also provides a 
competitive observational constraint on EoR. Based on the 2012 
Hubble Ultra Deep Field galaxies at z ~ 7-12, Ellis et al. 
(2013) show a continuous decline in the abundance of star¬ 
forming galaxies over 6 < z < 10, whereas Dunlop et al. (2013), 
based on the constancy of the UV continuum slope measured in 
z - 7-9 galaxies over a wide range in luminosity, show that stars 
at these redshifts are already enriched by earlier generations, 
supporting the idea of a reionisation process that extended be¬ 
yond z = 9. Finally, Robertson et al. (2013) combined the results 
from UDF12 with the measured Thomson optical depth and stel¬ 
lar mass density measurements (Stark et al. 2013) and provided 
constraints on the role of high-redshift star-forming galaxies in 
the reionisation process. Under reasonable assumptions about 
the escape fraction of ionising photons and the IGM clumping 
factor, they find that the observed galaxy population {Muy < 
-17 at z ~ 8) cannot simultaneously reionise the universe by 
z ~ 6 and produce a large Thomson optical depth, unless the 
abundance of star-forming galaxies and/or the ionising photon 
escape fraction increase in the range z ~ 15-25 with respect to 
what is currently observed (see e.g. Kuhlen & Faucher-Giguere 
2012). The star formation rate (SFR) history allows the dura¬ 
tion of reionisation to be probed (Madau & Dickinson 2014) be¬ 
cause star-forming galaxies provide most of the ionising pho- 



Fig. 1. Upper panel: Reionisation histories compared to low- 
redshift probes. Shaded envelopes show the low/high-r values 
derived from Planck with uncertainties. The red line indicates 
the best Planck value r = 0.66. Dashed lines are derived from the 
best ionisation rate (Vion) of ISHI15, ROB15, BOU15. Symbols 
are current constraints from direct line-of-sight measurements 
compiled by BOU15. Middle panel: Corresponding psfr for 
Planck and the best fit of ROB 15 with observational data com¬ 
piled by ROB 15. The grey dashed part corresponds to 2 hii ~ 1- 
Fower panel: Evolution of the integrated optical depth for the 
different studies plotted in the upper panel. 


tons (Robertson et al. 2013). Recent works by ISHI15, ROB 15, 
and BOU15 (see also references therein) have shown that star¬ 
forming galaxies at z ~ 8 - 10 were abundant enough to ionise 
the Universe and yield a Thomson optical depth of the order of 
0.06. 

CMB data can also constrain reionisation history through 
secondary anisotropies in temperature and polarisation. 
Measuring the latter has proven to be a challenging task because 
the signal is strongly contaminated by both foregrounds and 


2 

















Douspis et al.: A new parameterization of the reionisation history 


systematics on large angular scales (Planck Collaboration XXX 

2014) . This has affected the accuracy of optical depth mea¬ 
surements (Spergel et al. 2003; Page et al. 2007) and limited 
its precision (Hinshaw et al. 2013; Planck Collaboration XIII 

2015) . Recent results from Planck combining temperature and 
polarisation provide r = 0.078 + 0.019 in the A-CDM model 
(PCP15). The degeneracy between As and r can be broken by 
adding information about lensing (Planck Collaboration XVI 
2014). Planck also provides r = 0.66 + 0.013 when considering 
temperature, polarisation, lensing, and baryonic acoustic oscil¬ 
lations (BAO) in A-CDM (Eq. 17e of PCP15, hereafter referred 
to as the Planck reference). For an instantaneous reionisation, 
this corresponds to a redshift of Zre = ^-^- 12 ' Extending the 
standard A-CDM model probes reionisation history in a broader 
range of cosmologies. Two extreme cases for r are referred 
to in the following as high-r and low-r, corresponding to a 
MCMC run using both temperature and polarisation at high and 
low multipoles leaving the amplitude of lensing free and a run 
with temperature only at high multipoles and temperature and 
polarisation at low multipoles with curvature free, respectively ^ 
The two corresponding optical depth values are: r ~ 0.11 + 0.03 
and 0.054+ 0.021. 

We show in Fig. 1 (upper panel) the current data and con¬ 
straints from direct observations of the ionised hydrogen frac¬ 
tion 2 hii (symbols, from BOU15 and references therein), CMB 
(solid lines), and star-forming high-redshift galaxies (dashed 
lines); namely, we plot the measurements from line-of-sight ob¬ 
servations as compiled in ROB 15. We further add the lower limit 
estimates of the neutral hydrogen fraction at z ~ 7.7 from two 
additional Ly-a emitters (Faisst et al. 2014) and the neutral frac¬ 
tion estimate from the optical afterglow of the Swift gamma-ray 
burst (GRB) 140515A at z ~ 6.15 (Chornock et al. 2014). The 
high-T constraint from CMB (blue shaded area) is incompati¬ 
ble with other reionisation constraints. We furthermore note that 
within all the possible CMB models, those providing us with low 
optical depths agree best with the data points. However, changes 
in T only shift the best-fit CMB curves to high or low redshifts 
without modifying their shapes. As a result, no satisfactory fit 
seems attainable with the standard reionisation history used in 
CMB analyses. Figure 1 (middle panel) displays the evolution 
of the SFR density psfr compiled by ROB 15. Overplotted are 
the best models from Robertson et al. (2015) (dashed line) and 
the Psfr derived from the best-fit CMB models^ with increasing 
values of r from left (solid green curve) to right (solid light-blue 
curve). The dashed grey line shows the regime where the ionisa¬ 
tion fraction is set to unity, which is when ionisation is complete. 
Finally, in Fig. 1 (lower panel) we plot the Thomson optical 
depth derived from the ionising background due to star-forming 
high-redshift galaxies as a function of redshift (dashed lines). 
We also display the optical depth from the best-fit CMB mod¬ 
els (solid lines), together with the one-sigma measurement from 
PCP15 (shaded areas). As noted in BOU15 and ROB15, star¬ 
forming high-redshift galaxies provide enough ionising back¬ 
ground to obtain a Thomson optical depth of about 0.06, compat¬ 
ible (on the lower side) with Planck standard ACDM constraints 
(r ~ 0.066 + 0.013). Here again, the high-r value (blue) is defini¬ 
tively at odds with other astrophysical constraints. 


^ MCMC runs named base_omegak_plikHM_TTJowTEB and 
baseJilens_plikHM_TTTEEEJowEB are available on Planck Legacy 
Archive: http://pla.esac.esa.int/pla/ 

^ To convert from Nion to psfr, we assume \ogiMfesSon)) = 52.44 
as in ROB 15. 


3. New parameterization of the reionisation 
history 

Current astrophysical constraints seem to favour low values of 
the optical depth, regardless of the probe, with a sharp reioni¬ 
sation transition between redshifts ~ 7 and ~ 6. Reionisation 
histories derived from the abundance of star-forming galaxies 
suggest an asymmetric behaviour of Qhii with redshift. On the 
numerical side, recent simulations reproducing large volumes of 
the Universe allow for multiple populations of ionising sources, 
including early Population III stars and self-regulated UV emit¬ 
ting sources (e.g. Ahn et al. 2012; Park et al. 2013). These sim¬ 
ulations show that in the most general cases, the reionisation his¬ 
tory described in terms of Xe{z) shows a z-asymmetric behaviour 
(see e.g. Fig. 3 in Park et al. 2013). 

So far, CMB analyses in turn have assumed a redshift- 
symmetric instantaneous reionisation history (implemented e.g. 
in the publicly available code CAMB, Lewis et al. 2000). This 
widely used reionisation parameterization is a tanh-based ex¬ 
pression 2 hii(^) oc ^1 + tanh(^Et^^^^-^^)y2. The two key pa¬ 
rameters, Zre and Jz, measure the redshift at which the reaches 
half its maximum (typically 1.08, not including the second reion¬ 
isation of helium that occurs at z ~ 3.5) and the duration of 
reionisation, respectively. In this parameterisation Zre is related 
one-to-one to r for a fixed short duration. As shown in the upper 
panel of Fig. 1, the shape of the tanh-based expression does not 
provide a good fit to the data. The middle panel of Fig. 1 shows 
that the psfr evolution derived from CMB using a than-based 
parameterization does not agree with the data, or barely agrees 
on a rather narrow range of redshifts. 

The current astrophysical constraints on reionisation suggest 
that tanh-based histories have difficulty reproducing the data 
even in the low-r regime. We propose a new parameterisation 
that allows for a first period where reionisation can be slow and 
progressive, which is attributable to the softer ionising photons 
produced by the first stars and primordial dwarf galaxies, and a 
second period that can be faster leading, by z ~ 6, to the com¬ 
pletion of hydrogen and first helium reionisation by quasars, for 
example, with harder ionising photons. This parameterisation 
is motivated by the behaviour of gnu as deduced from simu¬ 
lations and observations. At the end of reionisation. Fan et al. 
(2006) show that the Gunn-Peterson optical depth evolves as 
tgp oc (1 + zT^'^ . The neutral hydrogen fraction in turn evolves 
as tgp/(1 + z)^^^, so we consider 1 - gnu oc (1 + z)^ at low red¬ 
shifts below 6. Fan et al. (2006) also note that at redshifts 6 and 
above, tgp grows much more quickly, possibly oc (1 + z)^^^ Such 
an accelerated phase is commonly observed in numerical simu¬ 
lations (Park et al. 2013; Fernandez et al. 2014). We assume that 
this phase is in fact the end of an exponential phase. 

The parameterisation we propose here thus follows these 
two asymptotic behaviours: a polynomial of order 3 below some 
pivot redshift (z;,) and an exponential one above: 

if z<Zp 1 - Qmi(z) ^ (1 + zf (1) 

if z>Zp 2hii(z) oc exp(-A(l + z)). (2) 

The parameterisation is described by three free parameters: the 
pivot redshift Zp, the gnu-value at the pivot, and the evolution 
rate in the exponential. Proportionality coefficients can be de¬ 
rived from these three parameters. Depending on the redshift and 
the gnu-value at the pivot, the expression is able to reproduce a 
large family of reionisation histories. 

A forthcoming paper (Ilic et al., in prep.) describes the pa¬ 
rameterization and its effects on different observations in more 
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Fig. 2. Upper panel: Reionisation histories compared to low- 
redshift probes. Same as Fig. 1 with the magenta curve display¬ 
ing the best-fit from the parameterization. Lower panel: Same as 
Fig. 1 with the magenta curve showing the associated derived 

PSFR- 


detail, in particular the CMB. We show in the upper panel of 
Fig. 2 the same data points and models as in Fig. 1. We over¬ 
plot our parameterisation of the reionisation history fitted to 
the measured values, i.e. excluding lower and upper limits. The 
corresponding parameter values are: Zp = 6.1, QmiiZp) = 
0.99986, A = 0.73. Figure 2 (upper panel) shows that the new 
parameterisation performs as well as or slightly better than con¬ 
straints from ROB 15 and BOU15 (dashed). It outperforms the 
tanh-based reionisation expressions used in CMB analyses (solid 
lines) while having the same value of r as the Planck low model 
and thus quasi-identical EE angular power spectra. However, the 
redshift of reionisation as defined as QuiiiZreio) = 0-5 is differ¬ 
ent in the two approaches: Zreio = 7.05 for our parameterisa¬ 
tion and Zreio = 7.65 for the tanh expression. Furthermore, the 
parameterisation reproduces the results from BOU15 in an eco¬ 
nomical way, with only three free parameters. Using this best fit, 
we derived the SFR density psfr from the parameterisation and 
compared it to the results and compilation from ROB 15 and to 
equivalent CMB-derived values using the tanh expression (Fig. 2 
lower panel). In contrast to the latter, the new parameterisation 
agrees perfectly with the data points and with the ROB 15 results 
in the redshift range probed by star-forming high-redshift galax¬ 
ies: It reproduces a slow beginning of reionisation that acceler¬ 
ates to reach an ionisation fraction of about 1, around the pivot 
point and a slowly varying ionisation state until the present. 


4. Conclusion 

Recent CMB observations favour lower values of the optical 
depth, in agreement with values derived from ionising back¬ 
ground estimates from star-forming galaxies and low-redshift 
line-of-sight probes. While this value is weakly sensitive to the 
modelling of the reionisation history in CMB studies, the latter 
has so far been assumed to be a simple step transition between 
neutral and fully reionised hydrogen. This naive parameterisa¬ 
tion does not render the behaviour of Qhii that is observed or 
derived from other probes. We propose here a new parameterisa¬ 
tion motivated by simulations and observations of the reionisa¬ 
tion fraction. This new parameterisation reproduces a large fam¬ 
ily of reionisation studies in an economical way for CMB studies 
and furthermore yields a good fit to the line-of-sight measure¬ 
ments of the ionised fraction and of the SFR densities. Moreover, 
it provides a more realistic value of the redshift of reionisation. 
We thus call for implementing this parameterisation in future 
CMB analyses for reionisation. Systematic studies of our param¬ 
eterisation, in particular in CMB studies, will be presented in a 
forthcoming paper. 
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